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Polycrystalline samples of LaNi2(Gei_a;Pa;)2 (x — 0, 0.25, 0.50, 0.75, 1) were synthesized and 
their properties investigated by x-ray diffraction (XRD) measurements at room temperature and 
by heat capacity Cp, magnetic susceptibihty x, and electrical resistivity p measurements versus 
temperature T from 1.8 to 350 K. Rietveld refinements of powder XRD patterns confirm that 
these compounds crystallize in the body-centered-tetragonal ThCr2Si2-type structure (space group 
lA/mmm) with composition-dependent lattice parameters that slightly deviate from Vegard's Law. 
The p{T) measurements showed a positive temperature coefficient for all samples from 1.8 K to 
300 K, indicating that all compositions in this system are metallic. The low-T Cp measurements 
yield a rather large Sommerfeld electronic specific heat coefficient 7 — 12.4(2) mJ/molK'^ for x — Q 
reflecting a large density of states at the Fermi energy that is comparable with the largest val- 
ues found for the ylFe2As2 class of materials with the same crystal structure. The 7 decreases 
approximately linearly with x to 7.4(1) mJ/molK^ for x = 1. The x measurements show nearly 
temperature-independent paramagnetic behavior across the entire range of compositions except for 
LaNi2Ge2, where a broad peak is observed at « 300 K from xiT) measurements up to 1000 K 
that may arise from short-range antiferromagnetic correlations in a quasi-two-dimensional magnetic 
system. High-accuracy Fade approximants representing the Debye lattice heat capacity and Bloch- 
Gruneisen electron-phonon resistivity functions versus T are presented and are used to analyze our 
experimental Cp(T) and p{T) data, respectively, for 1.8 K < T < 300 K. The T-dependences of 
p for all samples are well-described over this T range by the Bloch-Griineisen model, although the 
observed p(300 K) values are larger than calculated from this model. A significant T-dependence 
of the Debye temperature determined from the Cp(r) data was observed for each composition. No 
clear evidence for bulk superconductivity or any other long-range phase transition was found for 
any of the LaNi2(Gei_a;Pa;)2 compositions studied. 

PACS numbers: 74.70.Xa, 02.60.Ed, 75.20.En, 65.40.Ba 



I. INTRODUCTION 

The search for high temperature superconductors in- 
tensified after the discovery of superconductivity at 26 K 
in the compound LaFeAsOi_a;Fa;.^'^ Subsequent studies 
revealed even higher superconducting transition temper- 
atures (Tc) upon replacing La with smaller rare earth 
elements,^'* yielding a of 55 K for SmFeAsOi_a;Fa;.^ 
These compounds crystallize in the primitive tetrago- 
nal ZrCuSiAs (1 Ill-type) structure with space group 
PA/nmm.^ They have alternating FeAs and RO {R — 
rare earth element) layers stacked along the c-axis. The 
Fe atoms form a square lattice in the ah plane and are 
coordinated by As tetrahedra, where the coordinating As 
atoms lie in planes on either side of and equidistant from 
an Fe plane. The undoped parent compounds show cou- 
pled structural and antiferromagnetic (AF) spin density 
wave (SDW) transitions,^"^ which are both suppressed 
upon doping by partially substituting O by F. Such 
doping results in a nonintegral formal oxidation state 
for the Fe atoms. This suppression of the long-range 
ordering transitions appears necessary for the appear- 
ance of high-temperature superconductivity, as in the 
layered cuprate high-Tc superconductors. ^''''^'-'^^""'^^ The 
1111-type self-doped Ni-P analogue LaNiPO becomes su- 
perconducting at temperatures up to about 4.5 K.^*'^^ 

Subsequently, the parent compounds 4Fe2As2 [A — 



Ca, Sr, Ba, and Eu) were investigated.^ They crystal- 
lize in the body-centered-tetragonal ThCr2Si2 (122-type) 
structure with space group lA/mmm and contain the 
same type of FeAs layers as in the 1111-type compounds. 
In addition, they show SDW and structural transitions 
at high temperatures^^"^^ that are similar to those seen 
in the 1111-type compounds and these are similarly sup- 
pressed upon substituting on the A, Fe and/or As sites. ^ 
Superconductivity again appears at temperatures up to 
38 K as these long-range crystallographic and magnetic 
ordering transitions are suppressed. ^^"^^ However, the 
conventional electron-phonon interaction has been cal- 
culated to be insufficient to lead to the observed high 
Tc's, and strong AF fluctuations still occur in these com- 
pounds above Tc even after the long-range AF ordering 
is suppressed.^ The consensus is therefore that the su- 
perconductivity in these high-Tc compounds arises from 
an electronic/magnetic mechanism rather than from the 
conventional electron-phonon interaction.^ 

The above discoveries motivated studies of other com- 
pounds with the 122-type structure to search for new 
superconductors and to clarify the materials features 
necessary for high-Tc superconductivity in the ^Fe2As2- 
type compounds. For example, the semiconducting AF 
compound BaMn2As2 contains local Mn magnetic mo- 
ments with spin S = 5/2 and Neel temperature Tn = 
625 K.'^""'^'^ We recently doped this compound with K to 
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form a new series of metallic AF Bai_a;Ka;Mn2As2 com- 
pounds containing the Mn local magnetic moments, "^^ 
but no superconductivity has yet been observed in this 
system. ^'^''^^ This may be because the Tn was not suffi- 
ciently suppressed by the K-doping levels used. Studies of 
122-type compounds in which the Fe and As in ^Fe2As2 
are both completely replaced by other elements have also 
been carried out. For example, LaRu2P2 becomes super- 
conducting at Tc = 4.1 K."^^ SrPd2Gc2 was recently found 
to become superconducting with Tc = 3.0 K,'^'' with con- 
ventional electronic and superconducting properties. 

The Fc-based phosphides not containing magnetic 
rare earth elements such as CaFe2P2,^^ LaFe2P2,"'° 
SrFe2P2,^° and BaFe2P2 (Ref. 41) show Pauli paramag- 
netic behavior. The 122-type Co-based phosphides ex- 
hibit varying magnetic properties. SrCo2P2 does not 
order magnetically, although it has a large Pauli sus- 
ceptibility X ^ 2x10"'^ cm^/mol that has variously 
been reported to exhibit either a weak broad peak at 
~ 110 K attributed to "weak exchange interactions bc;- 
tween itinerant electrons" (Ref. 40), or a weak broad 
peak at ~ 20 K attributed to a "nearly ferromagnetic 
Fermi liquid" (Ref. 42). LaCo2P2 orders fcrroniagnct- 
ically at a Curie temperature Tq « 130 K,^'^''''^ and 
CaCo2P2 is reported to exhibit A-type antiferromag- 
netism at Tn = 113 K in which the Co spins align 
ferromagnetically within the basal plane and antiferro- 
magnetically along the c-axis.^* These differing magnetic 
properties of the Co-based phosphides are correlated with 
the formal oxidation state of the Co atoms, taking into 
account possible P-P bonding.^"* Compounds whore the 
Co atoms have a formal oxidation state of +2, ~ +1.5, 
and ;$ +1 show no magnetic order, ferromagnetic order 
and antiferromagnctic order, respectively.*^ None of the 
above 122-type Fe or Co phosphides were reported to 
become superconducting. 

Among Ni-containing 122-type compounds, supercon- 
ductivity has been reported with Tc = 0.62 K in 
SrNi2As2,*^ Tc = 0.70 K in the distorted structure of 
BaNi2As2,'*® Tc = 1.4 K in the orthorhombically dis- 
torted structure of SrNi2P2 (Ref. 47) and Tc = 3.0 K in 
BaNi2P2.*^ The Pauli paramagnet LaNi2P2 is reported 
not to become superconducting above 1.8 K.*^ There are 
confficting reports about the occurrence of superconduc- 
tivity in LaNi2Ge2 with either Tc = 0.69-0.8 K,^°'^^ or 
no superconductivity observed above 0.32 K.^^ 

Several studies have been reported on the normal state 
properties of LaNi2Go2. de Haas van Alphcn (dHvA) 
measurements at 0.5 K indicated moderate band ef- 
fective masses m* /m^ = 1.2 to 2.7, where rUe is the 
free electron mass.^-^ Electronic structure calculations 
were subsequently carried out by Yamagami using the 
all-electron relativistic linearized augmented plane wave 
method based on the density-functional theory in the 
local-density approximation.^^ The density of states at 
the Fermi energy E-p was found to be large, 'D{E-f) = 
5.38 states/(eVf.u.) for both spin directions, arising 
mainly from the Ni 3(i orbitals, where f.u. means formula 



unit. This T>{E-p) is comparable to the largest values re- 
ported for the FeAs-based 122-typc superconductors and 
parent compounds.^ Three bands wore fomid to cross E-p, 
with two Fermi surfaces that were hole- like (0.16 and 
1.11 holes/f.u.) and one that was electron-like (0.27 elec- 
trons/f.u.), and therefore with a net uncompensated car- 
rier charge density of 1.00 holes/f.u. Thus for the hypo- 
thetical compound ThNi2Ge2 one can assign formal ox- 
idation states Th"*"*, Ni+^ and Ge""*. Then substituting 
trivalent La for tetravalent Th yields a net charge carrier 
concentration of one hole per formula unit. The electron 
Fermi surface is a slightly corrugated cylinder along the c- 
axis centered at the X point of the Brillouin zone, indicat- 
ing quasi-two-dimensional character, similar to the elec- 
tron Fermi surface pockets in the FeAs-based 122-type 
compounds.^ On the other hand, the two hole Fermi sur- 
faces are three-dimensional and are centered at the X and 
Z (or M, depending on the definition^) points of the Bril- 
louin zone. These calculated Fermi surfaces were found 
to satisfactorily explain the results of the above dHvA 
measurements,^^ including the measured band masses. 
Prom a comparison of the calculated T>{Ef) with that 
obtained from experimental electronic specific heat data, 
Yamagami inferred that many-body enhancements of the 
theoretical band masses are small. 

Hall effect measurements on single crystals of 
LaNi2Ge2 arc consistent with the occurrence of multi- 
ple electron and hole Fermi surfaces, with the weakly T- 
dependent Hall coefficients given by a positive (hole-like) 
value i?H ^ +3 x 10^^*' m'YC for the applied magnetic 
field H parallel to the a-axis and a negative (electron- 
hke) value i?H ~ -2 x 10"!° m^/C for H parallel to the 
c-axis.^* The thermoelectric power obtained on a poly- 
crystalline sample of LaNi2Ge2 is negative. 

Herein we report our results on the mixed system 
LaNi2(Gei_a;P2;)2- For a; = or 1, alternating La 
and NiCe or NiP layers, respectively, are stacked along 
the c-axis. We wanted to investigate whether any new 
phonomena occur with Ge/P mixtures that do not oc- 
cur at the endpoint compositions, such as happens when 
the parent FeAs-based compounds are doped/substituted 
to form high-Tc superconductors. In addition, Yam- 
agami's electronic structure calculations for LaNi2Ge2 
discussed above^^ indicated some similarities to the elec- 
tronic structures of the FeAs-based 122-type compounds. 

We report structural studies using powder x-ray 
diffraction (XRD) measurements at room temperature, 
together with heat capacity Cp, magnetic susceptibility 
X, and electrical resistivity p measurements versus tem- 
perature T from 1.8 to 350 K for five compositions of 
LaNi2(Gei_a;Pa;)2 with < a; < 1. Our low-T limit 
of 1.8 K precluded checking for superconductivity with 
Tc < 1 K reported for LaNi2Ge2,^°'^^ but we did find ev- 
idence for the onset of superconductivity below ^ 2 K in 
two samples of LaNi2P2 from both p{T) and x{T) mea- 
surements. However, it is not clear from our measure- 
ments whether this onset arises from the onset of bulk 
superconductivity or is due to an impurity phase. 
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Also presented in this paper is the construction of Pade 
approximants''^ for the Debye and Bloch-Griineisen func- 
tions that describe the acoustic lattice vibration contri- 
bution to the heat capacity at constant volume Cy{T) 
of materials and the contribution to the p{T) of met- 
als from scattering of conduction electrons from acous- 
tic lattice vibrations, respectively. These Pade approxi- 
mants were created in order to easily fit our respective 
experimental data using the method of least-squares, but 
they are of course more generally applicable to fitting the 
corresponding data for other materials. The Debye and 
Bloch-Griineisen functions themselves cannot be easily 
used for nonlinear least-squares fits to experimental data 
because they contain integrals that must be evaluated 
numerically at the temperature of each data point for 
each iteration. Several numerical expressions represent- 
ing the Bloch-Griineisen^^"^^ or Debye^^ functions have 
appeared. However, they replace the integrals in these 
functions with infinite series, use very large numbers of 
terms, and/or use special functions. These approxima- 
tions are therefore not widely used for fitting experimen- 
tal data. One paper presented a method for approxi- 
mating the Debye function using the Einstein model. 
This method is also of little use for fitting because it uses 
a different equation for each temperature range and it 
becomes inaccurate at low temperatures. However, as 
we demonstrate, the Debye and Bloch-Griineisen func- 
tions can each be accurately approximated by a sim- 
ple Pade approximant over the entire T range. To our 
knowledge, there are no previously reported Pade ap- 
proximants for either of these two important functions. 
The T-dependences of p for all samples discussed here 
are well-described by the Bloch-Griineisen prediction, al- 
though the observed /9(300 K) values are larger than cal- 
culated. A significant T-dependence of the Debye tem- 
perature determined from the Cp(T') data was observed 
for each composition. 

The remainder of this paper is organized as follows. 
An overview of the experimental procedures and appara- 
tus used in this work is given in Sec. H. The construction 
of the Pade approximants for the Bloch-Griineisen and 
Debye fun(;tions is described in Sec. HI and Appendix A. 
The structural, thermal, magnetic, and electrical resis- 
tivity measurements of the LaNi2(Gei_a;Pa:)2 system and 
their analyses are presented in Sec. IV and Appendices B 
and C. A summary and our conclusions are given in 
Sec. V. 



II. EXPERIMENTAL DETAILS 

Polycrystalline samples of LaNi2(Gei_a:Pa;)2 {x = 0, 
0.25, 0.50, 0.75, 1) were prepared using the high purity 
elements Ni: 99.9+%, P: 99.999+%, and Ge: 99.9999+% 
from Alfa Aesar and La: 99.99% from Ames Laboratory 
Materials Preparation Center. Stoichiometric amounts of 
La, Ni, and Ge were first melted together using an arc fur- 
nace under high-purity Ar atmosphere. The arc-melted 



button was flipped and remelted five times to ensure 
homogeneity. Next, the samples (except for LaNi2Ge2 
which was prepared following the general procediircs out- 
lined in Ref. 64) were throughly ground and mixed with 
the necessary amount of P powder in a glove box un- 
der an atmosphere of ultra high purity He. The; pow- 
ders were cold-pressed into pellets and placed in 2 mL 
alumina crucibles. The arc- melted button of LaNi2Ge2 
was wrapped in Ta foil. The samples were then sealed 
in evacuated quartz tubes and fired at 990 °C for fa 6 
d. Samples containing phosphorus were first heated to 
400 °C to prereact the phosphorus. 

After the first firing, the phase purities of the sam- 
ples were checked using room temperature powder x- 
ray diffraction (XRD) with a Rigaku Geigerflex pow- 
der diffractometer and CuKq radiation. The x-ray pat- 
terns were analyzed for impurities using MDI Jade 7. If 
necessary, samples were thoroughly regroimd and repel- 
letized (except for LaNi2Ge2 which was just rewrapped 
in Ta foil) and either placed back in alumina crucibles 
or wrapped in Ta foil and resealed in evacuated quartz 
tubes. Samples were again fired at 990 °C for 5-6 d. The 
LaNi2P2 sample was arc-melted with additional La and P 
in order to achieve a single phase sample. This may have 
been necessary because the compound may not form with 
the exact 1:2:2 stoichiometry. After arc-melting, part of 
the sample was annealed for 60 h at 1000 °C. Throughout 
this paper, the annealed LaNi2P2 sample will be referred 
to as X = 1.00a and the as-cast sample as x = 1.00b 
where x is the composition of LaNi2(Gei_a;P2;)2. As seen 
in the XRD patterns and fits in Sec. IV A, all final sam- 
ples were single-phase except for two samples showing 
very small concentrations of impurities. 

Rietveld refinements of the XRD patterns were carried 
out using the FullProf package. Magnetization mea- 
surements versus applied magnetic field H and temper- 
ature T were carried out using a siiperconducting quan- 
tum interference device (SQUID) magnetometer (Quan- 
tum Design, Inc.). Gel caps were used as sample holders 
and their diamagnctic contribution was measured sepa- 
rately and corrected for in the data presented here. 

The Cp(T) and p(T) measurements were carried out 
using a Quantum Design Physical Property Measure- 
ment System (PPMS). Samples for heat capacity mea- 
surements had masses of 15-40 mg and were attached to 
the heat capacity puck with Apiezon N grease for ther- 
mal coupling to the platform. The p{T) measurements 
utilized a four-probe ac technique using the ac transport 
option on the PPMS. Rectangular samples were cut from 
the sintered pellets or arc-melted buttons using a jew- 
eler's saw. Platinum leads were attached to the samples 
using EPO-TEK PlOll silver epoxy. The sample was 
attached to the resistivity puck with GE 7031 varnish. 
Temperature-dependent p measurements were recorded 
on both cooling and heating to check for thermal hys- 
teresis. No significant hysteresis was observed for any of 
the samples. In addition, the vibrating sample magne- 
tometer (VSM) option on the PPMS was used to measure 
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the high-T magnetization of the LaNi2Ge2 sample up to 
1000 K. 



III. FADE APPROXIMANT FITS TO THE 
BLOCH-GRUNEISEN AND DEBYE FUNCTIONS 

A Pade approximant g{T) is a ratio of two polynomi- 
als. Here we write these polynomials as series in 1/T 
according to 



No 



T 



N2 



^0 + ^ 



(1) 



The first one, two or three and last one, two or three 
in each of the sets of coefficients Ni and Di in g{T) can 
be chosen to exactly reproduce both the low- and high-T 
limiting values and power law dependences in T and/or 
1 /T of the function it is approximating. This is a very 
important and powerful feature of the Pade approximant. 
Then the remaining terms in powers of 1 /T in the numer- 
ator and denominator have freely adjustable coefficients 
that are chosen to fit the intermediate temperature range 
of the function. A physically valid approximant requires 
that there are no poles of the approximant on the positive 
real T axis. 




A. Bloch-Griineisen Model 

The temperature-dependent electrical resistivity due 
to scattering of conduction electrons by acoustic lat- 
tice vibrations in monatomic metals is described by the 
Bloch-Griineisen (BG) model according to^^ 



p(T) = 47^(eR) 



T 



Or/T 



(e^ - 1)(1 - e-=^) 



where 



h 

72 



4n^/^aMfcBeR_ 



dx 
(2) 

(3) 



Gr is the Debye temperature determined from resistiv- 
ity measurements, Ti is Planck's constant divided by 27r, 
riccu is the number of conduction electrons per atom, 
M = (atomic weight)/A^A is the atomic mass, Nx is 
Avogadro's number, a = (volume/atom)^/'^, fce is Boltz- 
mann's constant, and e is the elementary charge. These 
variables map a monatomic metal with arbitrary crystal 
structure onto a simple-cubic lattice with one atom per 
unit cell of lattice parameter a. To calculate 7?.(9r) in 
units of i7cm, one sets the prefactor (?i/e^) in Eq. (3) to 
4108.24 Vt in SI units and calculates the quantities inside 
the square brackets in cgs units so that the quantity in 
square brackets has net units of cm. If one instead has 
a polyatomic solid, one can map the parameters of that 
solid onto those of the monatomic solid described by the 
Bloch-Griineisen model as explained in Sec. IV D below. 



FIG. 1: (Color online) (a) Normalized Bloch-Griineisen elec- 
trical resistivity in Eq. (5) (one in every five points used for 
fitting are plotted, open circles) and the Pade approximant 
fit (solid red curve), (b) Residuals (Pade approximant value 
minus Bloch-Griineisen formula value), and (c) Percent error 
(residual divided by value of the Bloch-Griineisen formula), 
all versus temperature T divided by the Debye temperature 
Or. 



In practice, one fits the T dependence of an experi- 
mental p{T) data set by the BG model using an inde- 
pendently adjustable prefactor p{Qr) instead of 47?.(6r) 
in Eq. (2), because accurately fitting both the magnitude 
and T dependence of a data set cannot usually be done 
using a single adjustable parameter 0r. One therefore 
normalizes Eq. (2) by p{T = Gr). When T = Or, the 
integral in Eq. (2) is 



dx = 0.2 366159, 



(e--l)(l-e-) 



yielding 



p(efl) = 0.9 464 635 7^(eR). 



(4) 



Equations (2)-(4) then yield the normalized T- 
dependence of the BG function (2) as 



^(^) - 4 226259 1^— 

p(eR) ~ ^-''^'^^Ur 

en/T 



(e^ - 1)(1 



-dx. 



(5) 
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TABLE I: Values of the coefficients in the Pade approximant 
in Eq. (8) that accurately fits the normalized Bloch-Griineisen 
function in Eq. (7). 
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This T-depcndcncc is only a function of tlie dimension- 
less normalized temperature T/Or. Therefore we define 
normalized p and T variables as 



Pn(T„) 



p(0r) 



(6) 



T = ^ 



and Eq. (5) becomes 



p„(T„) = 4.226 259 
ri/Tn 

X ' 



/ 

Jo 



(e^-l)(l-e- 



-dx. 



(7) 



A set of pn{Tn) data points calculated from Eq. (7) is 
plotted in Fig. 1(a). These values were then used as a 
set of "data" to fit by a Pade approximant as described 
next. 

In order to construct a Pade approximant function that 
accurately represents Pn{Tn) in Eq. (7), the power law 
Tn dependences of the latter function must be computed 
at high and low temperatures and the coefficients of the 
Pade approximant adjusted so that both of these limiting 
T dependences are exactly reproduced (to numerical pre- 
cision). Then the remaining coefficients are determined 
by fitting the approximant to a set of data obtained by 
evaluating the function over the entire relevant tempera- 
ture range. This procedure is described in Appendix A 1, 
which constrains the values of Di, D2, D3, and Dg- The 
resulting approximant is 



Pn(rn) 



El 



N2 
T\7 



N3 
T\7 



(8) 



where the fitted coefficients in the Pade approximant (8) 
are listed in Table I. The denominator in Eq. (8) was 
checked for zeros and all were found not to lie on the 
positive Tn axis. This insures that the approximant does 
not diverge at any (real positive) temperature. As seen in 



Fig. 1(b), the difference between the BG function and the 
Pade approximant is less than 1 x 10^"' at any T in the 
range < T/6r, < 2. By construction, the Pade approx- 
imant must asymptote to the exact BG T dependences 
at high- and low-T, respectively. Figure 1(c) shows the 
percent error of the fitted approximant. This error is 
largest at low T with a value of « 3% at T/0r « 0.1. 
This is acceptable considering the very small value of pn 
at such low Tn. For the most accurate fit of low-T ex- 
perimental data by a power law in T, one would directly 
fit experimental data by the power law rather than using 
the Pade approximant function. 

When fitting experimental p{T) data by the Pade ap- 
proximant Pn{Tn) in Eq. (8), one fits only the T depen- 
dence and not the magnitude of p{T) by the BG theory, 
because as noted above, one cannot in general obtain a 
good fit of both the magnitude and the T dependence 
of a measured p{T) data set using only the single fitting 
parameter 6r,. Thus, one fits an experimental p{T) data 
set by 



p{T)=po + pien)PniT/ejO 



(9) 



where po is the residual resistivity for T — )■ 0. The mean- 
ing of /9n(T/6R) is that one substitutes T/0r for Tn, 
according to Eqs. (6), in the Pade approximant function 
Pn(Tn) in Eq. (8). The three adjustable parameters po, 
p(0r) and are varied independently to obtain a fit 
to the data. Once a good fit is obtained and all three pa- 
rameters are determined, one can compare the measured 
value of p{T = Or) with the value predicted by the BG 
theory in Eq. (4). Often the agreement is not very good 
even for s- or sp-metals.®^ 



B. Debye Model 

The Debye model^*^ is widely used for fitting experi- 
mental heat capacity Cp (T) data taken at constant pres- 
sure p arising from acoustic lattice vibrations. It is some- 
times useful to fit a large T range of experimental Cp data 
using the Debye function. Here we describe the construc- 
tion of an accurate Pade approximant function of T that 
can easily be used in place of the Debye function (10) 
for least-squares fitting experimental Cp(T) data over an 
extended T range. It can also be conveniently used to cal- 
culate the T dependence of the Debye temperature from 
experimental lattice heat capacity data over an extended 
T range. 

The lattice heat capacity at constant volume V per 
mole of atoms in the Debye model is given by^'' 



Cv(T) = 9E 



3 

D/ Jo 



(e^ - 1)' 



■ dx, 



(10) 



where R is the molar gas constant and ©d is the De- 
bye temperature determined from heat capacity measure- 
ments. The Cv and T can be normalized to become di- 
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FIG. 2: (Color online) (a) Plot of normalized Debye function 
in Eq. (12) (One in every two points used for fitting are plot- 
ted) (open circles) and Fade approximant (red line), (b) Plot 
of residuals (Pade approximant minus Debye function), (c) 
Percent error (residual divided by value of the Debye func- 
tion). 



mensionless according to 



T = 



Equation (10) then becomes 



Cn(rn) = 9T„^ 



Cv(rn) 

R '■ 

T 

1/T„ 



(11) 



(e^ - !)■ 



■ dx. 



(12) 



The Cn was calculated for a representative set of val- 
ues using Eq. (12). The resulting data are plotted as 
open red circles in Fig. 2(a). 

As discussed in Appendix A 2, the Pade approximant 
that fits both the high- and low-T power law asymptotics 
of Cn(Tn) in Eq. (12) and has additional terms in pow- 
ers of 1/T in the numerator and denominator to fit the 
intermediate T range is 



c„(r„) 



N2 



t7' 



D2 

t7' 



t7' 



D8 



(13) 



where the coefhcients are given in Table II. The result- 
ing fit and error analyses are shown in Fig. 2. The Pade 



TABLE II: Values of the coefficients in the Pade approximant 
in Eq. (13) that accurately fits the normalized Debye function 
in Eq. (12). 



Coefficient 


Value 


No 


226.684 46 


Ni 


64.752 051 1 


N2 


17.285 710 5 


Na 


1.052 246 63 


Ni 


-0.035 843 776 1 


iVs 


0.027 925 482 7 






Di 


21.5840170 


D2 


9.539 977 83 


Da 


1.427 2431 




0.337 538 084 


Ds 


0.034 609 046 3 


De 


0.007 440 025 83 


Dv 


-0.000 210 411972 


Ds 


0.000 119 451046 



approximant does not deviate from the normalized De- 
bye function in Eq. (12) by more than 2 x 10^^ at any 
T as seen in Fig. 2(b). By construction, the deviation 
goes to zero at both low and high T. The percent error 
in Fig. 2(c) has its maximum magnitude of 0.3% at low 
T, and occurs because Cn{Tn — )■ 0) ^> and numerical 
precision becomes an issue there. For another example 
of the high accuracy and use of this Pade approximant, 
see Fig. 9 below, where the 8d versus T is calculated 
directly for each of our samples of LaNi2(Gei_a;P2;)2 us- 
ing the Debye function in Eq. (10) and compared with 
that found using the Pade approximant in Eq. (13); only 
small differences are found. 

To fit experimental Cp{T) data by the Pade approx- 
imant Cn(Tn) in Eq. (13), one fits both the magnitude 
and T dependence of Cp simultaneously using 



Cp(T) - nRCniT/Qu), 



(14) 



where n is the number of atoms per formula unit and Od 
is the only fitting parameter. Here one substitutes T/6d 
for Tn, according to Eqs. (11), in the Pade approximant 
function Cn(Tn) in Eq. (13). For a metal, one can add to 
Eq. (14) a linear specific heat term 7T giving 



Cp{T)^-fT + nRCn{T/eu). 



(15) 



The 7 is the Sommerfeld electronic specific heat coeffi- 
cient that can be experimentally determined from a prior 
separate fit to Cp(T) data at low T according to^^ 



T 



1 + 



(16) 



as in Fig. 6 below, where j5T'^ is the low-T limit in 
Eq. (A8) of the Debye heat capacity. In Eq. (15), the 
only fitting parameter is 0d- 
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FIG. 3: (Color online) (a) Unit cell lattice parameters a 
and c and (b) unit cell volume Vcoii versus composition x 
for LaNi2(Gei_a:Pa:)2. The error bars are smaller than the 
symbol size and the solid curves are guides to the eye. 



IV. EXPERIMENTAL RESULTS AND 
ANALYSES 



Structure and Chemical Composition 
Determinations 



The starting parameters for the Rietveld refinements of 
the powder XRD patterns were those previously reported 
for LaNiaPa (Refs. 68, 69, 70) and LaNiaGez (Refs. 53, 
64, 71, 72) that are presented in Table III. All samples in 
the LaNi2(Gei_a;Pa;)2 system were found to crystallize in 
the body-centered-tetragonal ThCraSia structure (space 
group lA/mmm) as previously reported for the compo- 
sitions X = and 1, and our refined values for the lat- 
tice parameters are in agreement with reported values, 
as shown in Table III. The refinements of the powder 
XRD patterns are shown in Figs. 15-18 in Appendix B 
and the crystal data are listed in Table III. All samples 
were also refined for site occupancy and no significant 
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FIG. 4: (Color online) Comparison of a section of the room 
temperature powder XRD pattern of LaNi2P2 before and af- 
ter annealing. Numbers above the peaks are their (hkl) Miller 
indices. 



deviations were found from the value of unity. However, 
the Rietveld fits were not very sensitive to changes in site 
occupation. 

The lattice parameters a and c and the unit cell volume 
Vccw are plotted versus composition x in Fig. 3. As the 
concentration of P increases, the lattice parameters and 
unit cell volume all decrease monotonically while, from 
Table III, the c/a ratio increases. The composition de- 
pendences of the quantities in Fig. 3 deviate slightly from 
the linearities expected from Vegard's Law. From Ta- 
ble III, the zp /Ge c-axis position parameter of the P / Ge 
site has a small overall increase with increasing P con- 
centration. 

An interesting effect was observed in the x-ray data for 
LaNiaPa. The as-cast sample after arc-melting showed 
broadening of the diffraction peaks with large c-axis con- 
tributions. After annealing the arc-melted sample for 
60 h at 1000 °C, those peaks became sharp and shifted 
to lower 29 angles reflecting an increased c-axis lattice pa- 
rameter. These effects are shown on an expanded scale 
in Fig. 4 for the (105) and (116) reflections. The c-axis 
peak broadening may arise from disorder in the interlayer 
stacking distances along the c-axis.'''^ 



B. Heat Capacity Measurements 

Plots of Cp(T) for our samples of LaNi2(Gei_2;Pa;)2 
from 1.8 to 300 K are shown in Fig. 5. The heat capac- 
ity at 300 K ranges from 118.3-123.0 J/molK for these 
samples. These values are approaching with increasing T 
the expected classical Dulong-Petit value Cy = ZnR = 
124.7 mJ/molK for the heat capacity due to acoustic 
lattice vibrations, where n is the number of atoms per 
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TABLE III: Crystallographic parameters of the body-centered-tetragonal LaNi2(Gei_a;Pi)2 system at room temperatme (space 
group M/mmm). Atomic coordinates are: La: (0, 0, 0), Ni: (0, 1/2, 1/4), P/Ge: (0, 0, Zp/Gc)- Listed are the lattice parameters 
a and c, the c/a ratio, the unit cell volume T4eii, and the z-coordinate of the P/Ge site zp/cc- The quality-of-fit parameters 
-Rp, Rv/p and are also listed. 



Compound 


a(A) 


c(A) 


c/a 


Kcii (A^) 


Zp/Gc 


Rp (%) 


-Rwp (%) 




Ref. 


LaNii.70(i)P2 


4.018(2) 


9.485(6) 


2.361 


153.1 


0.3716(2) 








68 


LaNi2P2 


4.010(1) 


9.604(2) 


2.395 


154.5 


0.3700(2) 








69 




4.007 


9.632 


2.404 


154.6 










70 


(as cast) 


4.027d(J) 


9.5073(9) 


2.3605(4) 


154.22(4) 


0.3592(7) 


15.7 


20.9 


4.24 


This work 


(annealed) 


4.0145(2) 


9.6471(6) 


2.4031(3) 


155.47(3) 


0.3681(6) 


13.2 


17.7 


3.01 


This work 


LaNi2(Po.75Geo.25 


2 4.0550(1) 


9.7289(3) 


2.3992(1) 


159.97(1) 


0.3685(3) 


11.5 


16.4 


9.05 


This work 


LaNi2(Po.5oGeo.50 


2 4.08132(7) 


9.7424(2) 


2.38707(9) 


162.281(9) 


0.3678(2) 


10.8 


UA 


3.80 


This work 


LaNi2(Po.25Geo.75 


2 4.1353(4) 


9.8012(9) 


2.3701(4) 


167.61(5) 


0.3668(2) 


9.40 


13.4 


6.35 


This work 


LaNi2 Ge2 


4.18586(4) 


9.9042(1) 


2.36610(6) 


173.535(6) 


0.3678(1) 


10.2 


13.2 


8.66 


This work 




4.1860(6) 


9.902(1) 


2.366 


173.51 


0.3667(2) 








71 




4.187 


9.918 














53 
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9.918(10) 


2.369 


173.8 
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72 




FIG. 5: (Color online) Heat capacity Cp versus temperature T 
for samples in the LaNi2(Gei-a;Pa;)2 system (open symbols). 
Fits of the data by Eq. (15), which is the sum of electronic 
and lattice contributions, are shown as solid curves with the 
respective color. For clarity, each plot is offset vertically by 
15 J/molK from the one below it. 



FIG. 6: (Color online) Heat capacity Cp divided by tempera- 
ture T versus (open symbols) and linear fits (solid curves of 
corresponding color) to the lowest-T data by Cp/T = j + fiT^. 
Table IV lists the temperature ranges of the fits and the values 
of 7 and /? obtained. 



formula unit (n = 5 for our compounds). 

To determine the Sommerfeld electronic specific heat 
coefRcient 7 and a low temperature value of the Debye 
temperature Od for each sample, the lowest temperature 
linear Cp{T)/T data for each sample, plotted in Fig. 6, 
were fitted by Eq. (16) and the values of 7 and f3 ob- 
tained. A value of 9d can be calculated from each value 
of /3 using^'' 



The values obtained for 7, /?, and Od for each sam- 
ple are listed in Table IV. A plot of 7 versus x in 
LaNi2(Gei_a;Pa;)2 is shown in Fig. 7, where a nearly lin- 
ear decrease in 7 with increasing x is seen. Using this 7 
values in Table IV, we then fitted the data from 1.8 to 
300 K in Fig. 5 by Eq. (15) and obtained the Qu fitting 
parameters listed in Table IV. The fits are shown as the 
solid curves in Fig. 5, which are seen to agree rather well 
with the respective data. However, small T-dependent 
deviations between the data and fit for each sample are 
seen, which we address next. 

Deviations of a Debye model fit from experimental lat- 
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TABLE IV: Values of 7 and [3 obtained from the low-T fits of the data in Fig. 6 by Eq. (16) are hsted together with the density 
of states at the Fermi energy T>{Ef) in units of states/(eV f.u.) for both spin directions calculated from 7 using Eq. (24). Also 
shown are the 0d values calculated from the low-T /3 values using Eq. (17) and from a global fit to all the lattice Cp{T) data 
from 1.8 to 300 K for each sample. Available values from the literature are also listed. 



Sample 


Low-T 
Fit Range 
(K) 


7 

(m,]/molK2) 


T>{Ey) 
(eVf.u.)-^ 


(mJ/mol K-*) 


Od 

Low-T fit 
(K) 


Od 

AU-T fit 
(K) 


Ref. 


LaNi2P2 (as-cast) 


1.81-5.34 


7.7(2) 




0.086(7) 


483(14) 


369(2) 


This work 


(annealed) 


1.81-7.31 


5.87(2) 


2.49 


0.126(2) 


426(3) 


365(3) 


This work 


LaNi2(Po.75Geo.25)2 


2.93-8.60 


7.4(1) 


3.13 


0.159(3) 


394(3) 


348(2) 


This work 


LaNi2(Po.5oGeo.5o)2 


1.82-7.13 


9.3(2) 


3.94 


0.194(7) 


369(5) 


326(2) 


This work 


LaNi2(Po.25Geo.75)2 


2.59-6.52 


11.27(6) 


4.78 


0.261(2) 


333.9(9) 


301(2) 


This work 


LaNi2Ge2 


2.93-6.41 


12.4(2) 


5.26 


0.371(9) 


297(2) 


287(2) 


This work 






14.5 




0.273 


328 




52 






12.7 (calc) 


5.38 








53 






13.5 (obs)" 










53 



"No experimental evidence or reference citation was given for this 
quoted observed value. 



12 - 



LaNi (Ge P ) 

2 l-x X 2 



O 



0.0 



0.2 



0.4 



0.6 



1.0 



FIG. 7: (Color online) Sommerfeld electronic specific heat 
coefficient 7 versus composition x in LaNi2(Gei_a:Pi:)2 from 
Table IV. For x = 1, the datum in Table IV for the annealed 
sample is plotted. The error bars are smaller than the data 
symbols. The line is a guide to the eye. 

tice heat capacity data are due to the following assump- 
tions and approximations of the model. 

• The system is assumed to be at constant volume as 
T changes, rather than at constant pressure which 
is the experimental condition. This deficiency can 
be corrected for if the T-dependent compressibility 
and thermal expansion coefficient are known for the 
compound of interest. 

• A quadratic density of phonon states versus energy 
is assumed, which terminates at the Debye energy 
/cb^d- For actual materials, this assumption can 
only be accurately applied at temperatures T <^ 
0D, which gives the Debye law [the second term 



on the right-hand side of Eq. (16)]. 

• Assumptions are made that the speed of acous- 
tic sound waves in a material is temperature- 
independent, is isotropic and is the same for lon- 
gitudinal and transverse acoustic sound waves. In 
general, these assumptions are too simplistic and 
the parameters are temperature dependent. 

• The Debye model only accounts for acoustic lattice 
vibrations and does not take into account optic lat- 
tice vibrations arising from opposing vibrations of 
atoms with different masses in the unit cell. The 
contribution of these to Cp (T) can be modeled by 
adding Einstein terms to the fit function. 

Because of these approximations and assumptions of 
the Debye model, the lattice heat capacity of a material 
is never precisely described by the Debye model over an 
extended temperature range such as from 2 K to 300 K. 
The most serious approximation in our T range is the 
second approximation. Within the Debye model 8d is 
independent of T . One can therefore parameterize the 
deviations of a fit from the data by allowing 0d to vary 
with 

The value of Od was calculated for each data point 
in Fig. 5, after subtracting the contribution from 7r ac- 
cording to Eq. (15), using the Debye function in Eq. (10). 
The resulting T dependences of Od are shown in Fig. 8, 
where Od is seen to vary nonmonotonically and by up to 
30% with increasing T . The plots have a similar shape 
to that for sodium iodide. The Od is expected to be 
constant below Odo/50 and above Odo/2, where Odo 
is the zero-temperature value of Od-^"* Our OD(r) data 
qualitatively agree with these expectations. 

Considerable scatter in the &y){T) data in Fig. 8 oc- 
curs at temperatures above 250 K and also below 7 K 
(although not as clearly visible in the figure). In these 
T ranges, Cp is becoming nearly independent of T, so in 
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FIG. 8: (Color online) Debye temperature Od versus temper- 
ature T (open symbols) obtained by solving Eq. (10) for each 
Cp(T) data point after subtracting the electronic contribu- 
tion 7T. The error bars plotted are calculated using Eq. (18). 
Solid curves are guides to the eye. 



these T ranges any error in the value of Cy is greatly 
amplified when 0d is calculated. The error bars on the 
values of QYi{T) therefore increase significantly in these 
T regions. We now consider such errors and for this dis- 
cussion ignore the difference between Cy and Cp. 

The scatter dOn in the derived 9d versus T depends 
on the statistical error dCy in Cy, 



dCy/deu' 



(18) 



This expression was used to obtain the error bars plotted 
in Fig. 8. The denominator dCy/dQu was calculated 
using the Debye function in Eq. (10). 

In order to more clearly see why the error in 0d sub- 
stantially increases above w 250 K, the high-T approxi- 
mation in Eq. (A7) can be used, yielding 



Cy ~ 3i? 



3i?eD 



(T » Od) 



20T2 

Taking the derivative with respect to 0d gives 

Inserting this result into Eq. (18) gives the approximation 
/ 10T2 



dOr 



V3i?eD 



dCy. 



(T » Od) 



This result shows that the error in 0d is proportional to 
at high T, which results in a dramatic increase in the 
scatter and error in Qu{T) at high T. 
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FIG. 9: (Color online) Difference AOd between the Oo cal- 
culated from the Debye function [Eq. (10)] plotted in Fig. 8 
and calculated from the Fade approximant [Eqs. (13) and 
(14)]. Solid curves are guides to the eye. For clarity, each 
plot is offset vertically upwards by 0.2 K from the one below 
it. Horizontal dotted lines are at AQu = for the data set 
with the corresponding color. 



Using the low-T approximation in Eq. (A8) and the 
same procedure as described above, the error in 0d is 



dOr 



dCy. 



(T < en) 



Therefore, similar to the situation at high-T, a small er- 
ror in Cp at low T is greatly amplified when calculating 

eD(T). 

To verify the applicability of the Fade approximant for 
the Debye function developed in Sec. IIIB, the 6d(T) 
was calculated for each data point in Fig. 5 using the 
Fade approximant in Eq. (13) instead of by directly us- 
ing the Debye function in Eq. (10). The electronic 7T 
contribution was again subtracted from the Cp(T) data 
first. The difference between these <du{T) values and 
those calculated using the Debye function in Eq. (13) is 
plotted versus T in Fig. 9. The values calculated from the 
Fade approximant do not deviate by more than 0.35 K 
from those calculated using the Debye formula. This er- 
ror is of order 0.1% of 6d, which is usually small com- 
pared to the error in 0d itself and is negligible compared 
to its T dependence. Therefore the Fade approximant 
provides a viable alternative for calculating 0d(T) that 
does not require evaluation of the integral in the Debye 
function (10) or the use of a look-up table for each data 
point. 
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FIG. 10: (Color online) Magnetic susceptibility x versus tem- 
perature T for the LaNi2(Gei_3;Pa;)2 system, (a) Measured 
X = M/H data (uncorrected), (b) Intrinsic x obtained af- 
ter correcting for both ferromagnetic (FM) and paramagnetic 
(PM) impurities. Sohd symbols of corresponding shape and 
color are values of x found from fitting M{H) isotherms. 



Magnetization and Magnetic Susceptibility 
Measurements 



Magnetization versus applied magnetic field M{H) 
isotherms were measured for the LaNi2(Gei_2;Pa;)2 sys- 
tem for H ~ 0-5.5 T and the results are plotted in 
Figs. 19-21 of Appendix C. The M versus T for the 
samples were also measured from 1.8 to 300 K at a fixed 
field H — iT and the resulting susceptibilities x = M/H 
are plotted in Fig. 10(a). As evident from the nonlinear 
behavior at low fields in the M{H) plots and the upturn 
that follows the Curie- Weiss-like behavor [x = C/ {T — 9)] 
in the observed x{T) at low temperatures, we infer the 
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FIG. 11: (Color online) Expanded plot along the vertical 
axis of the magnetic susceptibility x versus temperature T 
for LaNi2Ge2 up to 1000 K after correction for paramagnetic 
and ferromagnetic impurity contributions. The data below 
350 K were measured with a SQUID magnetometer and the 
data above 300 K were measured with a VSM, both at ap- 
plied fields of 3 T. Inset: Expanded plot of the SQUID data 
at low temperatures. The small downturn below about 10 K 
is most likely spurious, arising from an imperfect correction 
for the paramagnetic impurities. 



presence of saturating paramagnetic and/or ferromag- 
netic impurities in the samples. In order to determine 
the intrinsic behaviors, it is necessary to correct for these 
impurities. To determine the individual contributions to 
the susceptibility data, the M{H) curves were fitted by 



M(T, H) = Mo +xH + /M,atSs 



, (19) 



where Mq is the saturation magnetization of the ferro- 
magnetic impurities, x is the intrinsic susceptibility of 
the sample, / is the molar fraction of paramagnetic im- 
purities, g is the spectroscopic splitting factor of the im- 
purities which was fixed at g = 2 to reduce the number 
of fitting parameters, /zb is the Bohr magneton, fee is 
Boltzmann's constant, 9 is the Weiss temperature of the 
paramagnetic impurities (included for consistency with 
a possible Curie- Weiss law behavior at low H/T), M^at 
is the saturation magnetization of the paramagnetic im- 
purities, and Bs is the Brillouin function. The Brillouin 
function is 



Bs{x)^—{{2S+l)coth [(25+1) I 
where 

5MB 



coth 



(20) 
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FIG. 12: (Color online) Zero-field-cooled (ZFC) low- 
temperature measurement of the magnetic susceptibility x 
versus temperature T for LaNi2P2 with applied field H — 
20 G. 

and the molar saturation magnetization Mg^t of the para- 
magnetic impurities is 

where S is the spin of the impurities and A^a is Avo- 
gadro's number. 

To determine the values of /, 5*, and 9 for the param- 
agnetic impurities, a global two-dimensional surface fit of 
the M{H) data from 1-5.5 T taken at both 1.8 and 5 K 
was done using Eq. (19) with MATLAB's Surface Fitting 
Tool. For this H range, the ferromagnetic impurities are 
expected to be nearly saturated, as assumed by Eq. (19). 
The parameter values obtained are given in Table V. Fix- 
ing the variables /, S, and 6 at the respective values for 
each sample, each M{H) curve at higher temperatures 
was fitted by Eq. (19) in the range 1-5.5 T to obtain val- 
ues for Mo(T) and x(r). The M{H) fits are shown in 
Figs. 19-21 in Appendix C. A plot of il/g versus T is pre- 
sented in Fig. 22 in Appendix C and the fitted values of 
the intrinsic susceptibility x are plotted as solid symbols 
in Fig. 10(b). 

In order to correct the M{T) data at if = 3 T in 
Fig. 10(a) for the paramagnetic impurities, the above val- 
ues of /, 5, and 6 obtained by fitting the M{H) isotherms 
at 1.8 and 5 K were inserted into the last term of Eq. (19) 
to calculate their contributions versus T at 7? = 3 T. 
These contributions were then subtracted from the re- 
spective M [T] data that were already corrected for ferro- 
magnetic impurities to obtain the intrinsic susceptibility 
of the samples as plotted in Fig. 10(b). Previously re- 
ported values of x for LaNi2Ge2 are 24 x 10~^ cm'^/mol 
at 300 K (Ref. 50) and 28 x 10"^ cm^/mol at 296 K 
(Ref. 76) . The former value is essentially the same as our 
value 23.2 x 10~^ cm'^/mol at 300 K, as seen more clearly 



in Fig. 11 below where our x{T) data for LaNi2Ge2 are 
plotted on an expanded vertical scale for temperatures 
up to 1000 K. 

The x{T) data in Fig. 10 show that the samples in the 
LaNi2(Gei_2:Pa;)2 system exhibit nearly temperature- 
independent paramagnetism over the composition region 
X = 0.25-1. This trend does not extend to LaNi2Ge2 
(x = 0), for which a broad maximum appears to occur 
at K, 300 K, which is close to the upper temperature 
limit of the SQUID magnetometer. In order to deter- 
mine whether a maximum in x near 300 K does occur, 
M{H, T) data for the sample of LaNi2Ge2 were measured 
using a VSM from T = 300 to 1000 K. The M{H) 
isotherm data are plotted in Fig. 21(b) of Appendix C 
and the x{T) = M{T)/H data at iJ = 3 T are plotted in 
Fig. 11. The magnetic contributions due to the sample 
holder and paramagnetic impurities are corrected for in 
the plots. 

The data in Fig. 11 clearly show that a broad peak 
occurs in x{T) of LaNi2Ge2 at about 300 K. The plot in 
Fig. 10(b) shows the peak plotted on a less expanded ver- 
tical scale. The cause of this peak is not clear, but may be 
due to low-dimensional antiferromagnetic correlations.^^ 
Further investigation is needed. The inset in Fig. 11 
shows an expanded view of the low temperature behav- 
ior. Due to its smooth nature, the small downturn in the 
data below about 10 K is most likely spurious due to a 
slight error in correcting for the susceptibility contribu- 
tion of the paramagnetic impurities in the sample. 

The onset of superconductivity was observed in both 
the annealed and as-cast samples of LaNi2P2 at 2.6 and 
2.2 K, respectively, from zero-field-cooled (ZFC) magne- 
tization measurements in a field of 20 Oe. Figure 12 
shows these low-field x{T) = M{T)/H measurements. It 
is not clear whether the data represent the onset of bulk 
superconductivity in LaNi2P2 or if the diamagnetism 
arises from a superconducting impurity phase. The p{T) 
data in Fig. 13 below do not clarify this issue. 

We now analyze the normal-state x of the samples. 
The magnetic susceptibility of a metal consists of the 
sum of the spin and orbital contributions 

^^^spin^^orb^ (21) 

In the absence of local magnetic moments, the spin con- 
tribution x'^P™ is the Pauli susceptibility x^*^"'' of the 
conduction electrons. One can estimate ^p^,"'' using^ 

^Pauli ^ ^^2p(^p)^ (r = 0) (22) 

where g is the spectroscopic splitting factor and T>{E-p) 
is the density of states at the Fermi energy E-p. Setting 
g = 2 gives 

^Pauii ^ (3 X 10-^)I?(£;f), (23) 

where t^p^"'' is in units of cni"^/mol, 'D{Ep) is in units 
of states/eVf.u. for both spin directions and f.u. means 
formula unit. 
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TABLE V: Values obtained by simultaneous fitting of the 1.8 and 5 K M{H) isotherms. The parameters listed are / [molar 
fraction of paramagnetic (PM) impurities], S (spin quantum number of the PM impurities), 6 (Weiss temperature of the PM 
impurities), Afo (saturation magnetization of the ferromagnetic impurities), and \ (intrinsic magnetic susceptibility of the 

compouud). Tli(> negative signs of iiidicatc aiitifV-rrcjiuagiictic iiitt^ractitnis l)('tw(x>ii iiiagiitddc impurities. 





/ 


S 


e 


Mo 


X 


Sample 


(io-«) 




(K) 


(GcmVmol) 


(lO"'^ cm^mol) 


LaNi2P2 (annealed) 


3.89 


3.48 


-5.88 


0.0478 


2.92 


LaNi2 (Po. 75 Geo.26 ) 2 


2.33 


1.92 


-1.02 


0.1175 


5.53 


LaNi2(Po.BoGeo.6o)2 


4.12 


2.42 


-1.89 


0.4128 


7.49 


LaNi2(Po.25Geo.76)2 


5.73 


2.57 


-1.54 


0.1834 


16.0 


LaNi2Ge2 


2.35 


2.43 


-1.08 


0.0700 


19.5 



TABLE VI: Contributions to the magnetic susceptibility X- The parameters listed are the observed value x°'"' of X averaged 
over the temperature range measured and the contributions from the Pauli spin susceptibility x''*'"'' and the orbital core 
susceptibility x'^"'^'^) Landau susceptibility T^^andau^ Van Vleck susceptibility x^^- AH values are in units of 10~^ cm^/mol. 



Sample 




T^Pauli 


^core 


^Landau 


X^v 


LaNi2P2 (as-cast) 


1.49 


10.5(3) 


-18.8 


-3.50(1) 


13.3(4) 


LaNi2P2 (annealed) 


2.82 


8.04(3) 


-18.8 


-2.68(1) 


16.3(4) 


LaNi2(Po.75Geo.25)2 


5.28 


10.1(2) 


-19.3 


-3.37(7) 


17.9(3) 


LaNi2(Po.5o Geo. 50)2 


7.32 


12.7(3) 


-19.7 


-4.2(1) 


18.5(4) 


LaNi2(Po.26Geo.75)2 


14.36 


15.44(8) 


-20.1 


-5.15(3) 


24.7(1) 


LaNi2Ge2 


21.01 


17.0(3) 


-20.5 


-5.7(1) 


30.4(4) 



To calculate x^'""^'', one can obtain an estimate of 
T>{Ef) from the Somm(;rfc!l(i electronic linear specific 
heat coeSicient 7 according to 



where®^ 



70 



7 = 7o(l + Aep), 



^-V{Ef) = 2.359 V{Ef) 



(24) 



is the bare Sommerfeld coefficient in the absence of 
electron-phonon coupling, Aep is the electron-phonon 
coupling constant, and in the right-hand equality of 
Eq. (24) D{E-p) is in units of states/eVf.u. for both 
spin directions and 70 is in units of mJ/molK^. Tak- 
ing Aep = yields 

^Pauii ^ (^ 370 ^ 10-5)7, (25) 

where y^^^^^ is in units of cm^/mol and 7 is in units of 
mJ/molK2. The values of V{Ey) and x^''"" obtained 
using Eqs. (24) and (25), respectively, and using the 7 
values given in Table IV, are listed in Table VI. 

The orbital susceptibility consists of three contribu- 
tions 



.orb 



VV 



.Landau 



(26) 



where x'^°'^*' is the diamagnetic contribution from the 
atomic core electrons, is the paramagnetic Van Vleck 
susceptibility, and ^^Landau ^.j^g diamagnetic Landau 
susceptibility of the conduction electrons. The Landau 
susceptibility is approximated by^^ 



X 



Landau 



1 /TOe\ 

3 Vm*/ 



X 



Pauli 



For our samples, it is assumed that m* = mo, therefore 
.^Landau ^..^-^^ Calculated from the x^**"'' obtained above. 
An estimate of x'^"'^* was obtained from the sum of the 
Hartree-Fock diamagnetic atomic susceptibilities.^^ Us- 
ing the observed average value of the susceptibility (x"^*^) 
over the temperature range 2-300 K and the calculated 
sum of the values x^^"", X™'*' and x'^^'^'^'^", the Van Vleck 
susceptibility is calculated from Eqs. (21) and (26) via 



X 



vv . 



X 



obs 



Pauli 



_j_ •1^^°''^ _|_ T^Landau^ 



Pauli 



The values obtained in the manner described for x 

^core^ ^VV ^Landau yisied in Table VI. 



D. Electrical Resistivity Measurements 

The p versus T data for our LaNi2(Gei_2;P2;)2 sam- 
ples are plotted in Fig. 13(a). The magnitudes of the 
data may not be reliable due to the polycrystalline na- 
ture of the samples and the resulting grain boundary 
scattering. Such scattering is expected and found to be 
smallest for the samples of LaNi2Ge2 and LaNi2P2 that 
were cut from arc-melted buttons and therefore had a 
higher density than the other samples cut from sintered 
pellets. The values of p at 1.8 and 300 K are listed in 
Table VII along with previously reported values from the 
literature.^^'^^'^°'^^ The positive slopes of the p{T) data 
indicate that all samples are metallic. The maximum re- 
sistivities at 1.8 K and at 300 K both occur for x = 0.50 
where the disorder on the Ge/P sublattice is largest, as 
might have been anticipated, so this may be an intrinsic 
effect. 
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TABLE VII: Values of the Debye temperature determined from resistivity p measurements (6r), the resistivity at the Debye 
temperature [p(6r)] and the residual resistivity (po) obtained from least-squares fits of the p(T) data in Fig. 13(a) by Eq. (9). 
Also listed are p values at ~ 2 K and at 300 K, and literature values parallel to the c-axis (pc) and parallel to the o-axis {pa) 
of a single crystal. The systematic errors in our p values due to uncertainties in the geometric factors are of order 10%. 



Sample 


Or (K) 


p(eR) 


Po 


p(~ 2 K) 


p(300 K) 


Ref. 






(/iSl cm) 


(/jS7 cm) 


{^lQ. cm) 


(pO cm) 






91 1 ('9^ 




OO. lO^O j 


OO 


1 AR 


J_ lilb WUIK. 


(annealed) 


265(3) 


109.(1) 


26.0(1) 


25 


152 


This work 


LaNi2(Po.75Geo.25)2 


242(1) 


85.8(4) 


191.51(5) 


191 


300 


This work 


LaNi2(Po.5oGeo.5o)2 


208(1) 


102.5(7) 


249.79(9) 


249 


401 


This work 


LaNi2(Po.25Geo.75)2 


119(7) 


37.(2) 


95.9(3) 


95 


189 


This work 


LaNi2Ge2 


148(5) 


38.(1) 


6.8(2) 


6.1 


85 


This work 










0.4 




51 










~ 1 


~ 26 (p„), ~ 38 (pc) 


80 










~ 5 


~ 80 


81 










~ 1-2 


~ 43 


55 



As seen in the expanded plot at low temperatures down 
to our low-T limit of 1.8 K in Fig. 13(b), there appears 
to be an onset of superconductivity occurring at « 2.1 K 
for both the annealed and as-cast samples of LaNi2P2, 
consistent with the above measurements of the low-field 
magnetic shielding susceptibilities in Fig. 12. It is not 
known if these results indicate the onset of bulk super- 
conductivity or whether they are due to a superconduct- 
ing impurity phase. 

The T dependence of p was fitted by Eq. (9), which 
includes our Padc approximant for the Bloch-Griincisen 
function. From these fits shown in Fig. 13(a), estimates of 
the Debye temperature (Or), p at the Debye temperature 
[p(Qr)], and residual resistivity (po) were obtained as 
listed in Table VII. The Gr for the LaNi2(Po.25Geo.75)2 
sample is lower than the values for the adjacent compo- 
sitions. We speculate that this arises from inaccuracies 
introduced by the polycrystalline nature of the samples. 

The 6r obtained from p(T) measurements is usually 
different from that obtained by heat capacity measure- 
ments (9d) although in some cases agreement found. ''^ 
For oiu" samples in the LaNi2(Gci_2:PrE)2 system, while 
the magnitudes of 0r differ considerably from sample 
to sample, both sets of 9r and 6d data in Tables VII 
and IV, respectively, indicate an overall decrease in the 
Debye temperature with increasing Ge content. 

Like the Debye model, the BG model makes several 
assumptions and approximations in order to obtain an 
analytic formula for p{T): 

• A strong approximation is that the lattice vibra- 
tions have a Debye spectrum. 

• Another strong approximation is that Umklapp 
conduction electron scattering processes are ig- 
nored. Umklapp scattering is expected to be T- 
dependent in our T range. ®^ 

• Only longitudinal phonons arc assumed to con- 
tribute to electron-phonon scattering. 



• The calculation is carried out at constant volume 
instead of at constant pressure as in experiments. 

• The conduction electron energy is assumed not to 
change due to scattering off phonons. 

• The phonons are assumed to be in thermal equilib- 
rium. 

As a result of these assumptions and approximations, the 
T dependence of 0r is much stronger than that of 0d-^^ 
Zimon concluded, "The actual observed value of 0r does 
not have any great significance." Thus 0r is less mean- 
ingful than 6d- Furthermore, differences between Or 
and Gd are not unexpected, given the different assump- 
tions of the BG and Debye models enumerated here and 
in Sec. IV B, respectively. 

Now we compare our experimental result for the mag- 
nitude of p(300 K) for the annealed sample of LaNi2Ge2 
with the value predicted by the Bloch-Griineisen theory. 
Wc chose to use the p data for this compoimd for com- 
parison with the theory because it has the lowest resid- 
ual resistivity of our five samples. The BG theory de- 
scribes monatomic materials. In order to use this theory 
to predict the resistivity of polyatomic compounds such 
as LaNi2Ge2, it is necessary to slightly modify Eq. (3) as 



^ 1 4"-ccnafcB0R 



(27) 



where the variables have the same meaning as before, 
except 



iVA 

n 



n ^ 



i=l 



Mi 



(28) 



is the average inverse mass of the atoms in a f.u., n is the 
number of atoms per f.u. and is the atomic weight of 

element i. In this case, a = {Vcdi/nZy^^ because there 
are Z = 2 f.u. of LaNi2Ge2 with n = 5 per body-centered 
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FIG. 13: (Color online) (a) Electrical resistivity p versus tem- 
perature T for the LaNi2(Gei_2,P3;)2 system (open symbols). 
Solid curves of corresponding color are least-squares fits by 
Eq. (9), which includes the Pade approximant in Eq. (8) that 
is used in place of the normalized Bloch-Griineisen function 
in Eq. (7). Values of the fitting parameters obtained are listed 
in Table VII. (b) Expanded plots at low temperatures of p{T) 
for as-cast and annealed LaNi2P2. 



tetragonal unit cell with volume T^oii = o^c where a and 
c are given in Table 111. 

The number of carriers per f.u. in LaNi2Ge2 predicted 
by band structure calculations^'^ is 1.54. Therefore, 
nceW = 1.54/5 = 0.308 carriers/atom because the car- 
riers are modeled as being evenly distributed among the 
5 atoms in the f.u. As described previously, to calculate 
7?-(Br) in units of ficm in Eq. (27), the prefactor Ti/e^ 
is set to 4108.24 i7 and the part in the square brackets 
is calculated in cgs units. Using 0r — 148 K from Ta- 
ble Vll, 7^(8r) is calculated to be 3.763 ^.Q. cm. Inserting 
this result into Eq. (2) for T = 300 K, p(300 K) is pre- 
dicted to be 7.53 fiVlcm. This value is about a factor of 



10 smaller than the experimental value of 79 /ifi cm [after 
subtracting p(1.8 K)] from Table Vll, which is a large dis- 
agreement in magnitude even though the T dependence 
for p for LaNi2Ge2 was well-fitted by the BG theory. This 
large disagreement is typical of the BG theory as applied 
to transition metals and alkaline-earth metals^® and illus- 
trates why it is necessary to use an adjustable prefactor 
p(0r) in Eq. (9) in addition to the fitting parameter 0r 
in order to fit both the magnitude and T dependence of 
experimental data using the Bloch-Griineisen model. 

The reason for the enhanced electron-lattice resistivity 
in transition metals with both s- and d-bands crossing 
the Fermi energy is that the s-electrons carry most of 
the current due to their much lower effective mass than 
the d-electrons, and electron-phonon scattering from an 
s-band into a much higher density of states d-band acts 
as a conduction electron sink, thus strongly enhancing 
the resistivity due to electron-phonon scattering. 



V. SUMMARY AND CONCLUSIONS 

Single- or nearly single-phase polycrystalline samples 
of LaNi2(Gei_a;Pa;)2 with the compositions a; = 0, 0.25, 
0.50, 0.75, and 1 were synthesized and their properties 
measured. Rietveld refinements of the powder XRD 
patterns showed that all samples have the tetragonal 
ThCr2Si2 structure with space group lA/mmm. The 
refined crystal data are presented in Table III. A possi- 
ble stacking disorder was observed along the c-axis in the 
as-cast sample of LaNi2P2 as revealed by the c-axis line 
broadening in Fig. 4. 

Electrical resistivity p measurements showed a posi- 
tive temperature coefficient for all samples from 2 K 
to 300 K, indicating that all compositions in this sys- 
tem are metallic. Consistent with this result, the low-T 
Cp measurements yield a rather large Sommerfeld elec- 
tronic specific heat coefficient 7 = 12.4(2) mJ/molK^ 
for a; = refiecting a large density of states at the 
Fermi energy comparable to the largest values found for 
the ylFe2As2 class of materials with the same crystal 
structure.^ The 7 decreases approximately linearly with 
X to 7.4(1) mJ/molK^ for a; = 1. 

New Fade approximants for the Debye and Bloch- 
Griineisen functions were presented. They have the dis- 
tinct advantage of being able to precisely reproduce the 
power law dependences at both the low- and high-T lim- 
its of the function they are approximating. The Fade 
approximants presented here for each of Cy{T) and p(T) 
cover the entire T range and have a good balance of high 
accuracy and low number of terms. The T dependences of 
p for all samples were well-fitted by the Bloch-Griineisen 
model and values of the Debye temperature Or were ob- 
tained (Table VII), although the measured magnitudes 
were larger than calculated on the basis of this model. 
Fitting the T dependences of Cp revealed significant T 
dependences of the Debye temperatures 0d- 

Magnetization and magnetic susceptibility measure- 
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ments of our LaNi2(Gei_2:Pa;)2 samples revealed nearly 
temperature-independent paramagnetic behavior (except 
for LaNi2Ge2), with increasing susceptibility as the Ge 
concentration is increased. For LaNi2Ge2, a broad peak 
was observed in x(T) at « 300 K. A possible explanation 
is that the peak arises from the onset of strong antifer- 
romagnetic correlations in a quasi-two-dimensional mag- 
netic system/'^ There may be a correlation here between 
the x(^) behavior and the corrugated electron cylinder 
Fermi surface found by Yamagami^^ that was mentioned 
in the Introduction. However, the heat capacity data 
were well-fitted from 1.8 to 300 K for all samples by an 
electronic 7T term plus a lattice contribution described 
by the Debye model as shown in Fig. 5 and Table IV, with 
no clear evidence for an additional magnetic contribution 
in LaNi2Ge2. 

An interesting finding is the onset of a superconduct- 
ing transition at T « 2.2 K in both the annealed and 
as-cast samples of LaNi2P2. This onset was observed in 
both the resistivity and magnetic susceptibility measure- 
ments. However, it is not clear if these onsets are due to 
bulk superconductivity or to an impurity phase. Apart 
from the broad peak in x{T) for LaNi2Ge2 and the po- 
tential bulk superconductivity in LaNi2P2, there were no 
other signs of structural, magnetic, or superconducting 
transitions down to 1.8 K in the samples. 
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Appendix A: Pade Approximants 

1. Bloch-Griineisen Model 

In order to construct a Fade approximant function that 
accurately represents another function, the power law T 
dependences of the latter function must be computed at 
high and low temperatures and the coefficients of the 
Fade approximant adjusted so that both of these limit- 
ing T dependences are exactly reproduced (to numerical 
precision). For the EG function at high T, Qr/T <C 1 
in Eq. (7). Therefore, the integrand in Eq. (7) can be 
expanded in a Taylor series about .x = as 



(e^ - 1)(1 - e-^) 
Equation (7) then bcicomes 
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0{x''). {x <€. 1) 



p„(r„) « ATl+QT^ + BT^-^ +0{T-^) (T » Gr) 
A = 1.056 565 (Al) 
B = -0.058 698 04 



At low temperatures, Or/T — ^ og and the upper limit to 
the integral in Eq. (7) can be set to 00. The integral in 
Eq. (7) is then 



-dx = 5!C(5) « 124.431331, 



p x^ 
io (e--l)(l-e-)^ 

where C(z) is the Riemann zeta function. Inserting this 
result into Eq. (7) yields 

Pn = c^„^ (r<eR) (A2) 

C = 525.8790 

There are no additional terms in powers of Tn in the Tay- 
lor series expansion of pn{Tn) about Tn = 0, because it 
is not possible to express the exponentials in the integral 
in Eq. (7) as Taylor series in 1/x about 1/a; = 0. 

Based on the above low- and high-T expansions of the 
normalized EG function, the coefficients of the Fade ap- 
proximant are now chosen so that the limiting T depen- 
dences of the approximant exactly match the required 
power law T dependences in Eqs. (Al) and (A2). In 
addition, intcrniediate power terms were added in pairs 
(one in the numerator and one in the denominator) until 
there were enough for the approximant to accurately fit 
the intermediate temperature range of the EG function. 
The resulting approximant has the form 



Pn(Tn) = 



El. 



N2 



N3 



D2 



(8) 



In order that Taylor series expansions of the Fade ap- 
proximant (8) at high- and low-T correctly reproduce 
the coefficients of the limiting T dependences of the nor- 
malized BG function in Eqs. (Al) and (A2), some of the 
coefficients Ni and Di in Eq. (8) are not independent. Ex- 
panding Eq. (8) as a Taylor series in 1/Tn about l/Tn = 
gives 

, D2^No - D1D3N0 - D1D2N1 + Dl^N2 ^ _i 



(A3) 



Equating the coefficients in Eq. (A3) with the respective 
coefficients in the high-T Taylor series expansion of the 
normalized BG function in Eq. (Al) yields 

„ _ ^0 J. AN2-BN0 

The Taylor series expansion of the approximant (8) about 
Tn = is 



Pn = ^Tn5+C?(Tn*^). 



(A5) 



Equating the coefficient of this T^ term with that of the 
Bloch-Griineisen function in Eq. (A2) yields 



A3 
C 



(A6) 
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Using the four constraints in Eqs. (A4) and (A6), one can 
exactly reproduce the r„^, Tn^, and T^"^ dependences of 
the higli-T expansion in Eq. (Al) and the tJ' dependence 
of the low-T expansion in Eq. (A2) for the normaUzed BG 
function. 

A table of values generated from Eq. (7) and plotted 
above in Fig. 1(a) was least-squares fitted by Eq. (8) 
using the constraints in Eqs. (A4) and (A6) for Di, D2, 
D3, and Dg. The final values of the coefficients in the 
Pade approximant (8) are listed in Table I. 



2. Debye Model 

At high temperatures l/Tn ^ 1, the integrand in 
Eq. (12) can be expanded in a Taylor series about a; = 0, 
yielding 

4 r 4 

After evaluating the integral in Eq. (12) using this ap- 
proximation, the Cn(T„) becomes 

C7„(T„) = + OT-1 + FT„-2 + 0(T„-^) (T > Od) 

(A7) 

where 



and the corresponding Taylor series expansions of the 
Pade approximant (13) yield the constraints 



Do 



No 
E 

Ni 
E 

-FNo + EN2 



(A9) 



E-^ 



G ■ 



These constraints on Dq, Di, D2 and Dg, guarantee that 
the Tn°, Tn~\ and Tn ^ terms in the high temperature 
series expansion and the Tn'^ term in the low-T expansion 
of the Pade approximant (13) exactly match to within 
numerical accuracy the corresponding terms in the Tay- 
lor series expansion of the Debye function in Eqs. (A7) 
and (AS), respectively. 

The Debye fimction data in Fig. 2(a) were least-squares 
fitted by the Pade approximant in Eq. (13) using the con- 
straints in Eqs. (A9). After fitting, the denominator of 
the approximant was checked for zeros and all were found 
to be at non-real Tn. This insures that the approximant 
does not diverge at any (real positive) temperature. The 
final coefficients in the Pade approximant are listed in 
Table II. 



E = 3, 

In the limit of low temperatures l/Tn — >■ 00, the integral 
in Eq. (12) can be evaluated with an upper limit of 00 to 
obtain 



Appendix B: Rietveld Refinement Figures 

In this Appendix the Rietveld refinement Figs. 14-18 
referred to in Sec. IV A are shown. 



Cn(Tn) = 

G = 



GTn^ 



(T < Gd) 



(A8) 



This dependence of the lattice heat capacity at low 
temperatures is universal and is known as the Debye T'^ 
law. Similar to the Bloch-Griineisen function, it is not 
possible to obtain additional terms in the Taylor series 
expansion of Cn(Tn) in Eq. (12) about Tn = 0. 

Since the Pade approximant must follow Eqs. (A7) and 
(AS) at high and low temperatures, respectively, the ap- 
proximant was set up in the form 



Cn(Tn) = 



±n J n 



+ 



J- n 



+ ■ 



, Dj , Da 



(13) 



The number of intermediate power terms was increased in 

pairs (one in the numerator and one in the denominator) 
until the final fitted approximant accurately matched the 
Debye function in the intermediate T range. Using the 
same procedure as described in Sec. Ill A, the high- and 
low-T limits of the Debye function in Eqs. (A7) and (A8) 
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FIG. 14: (Color online) Room temperature powder XRD pat- 
tern (red open circles) of LaNi2P2, Rietfeld refinement fit 
(solid black line), difference profile (lower solid blue line), and 
positions of Bragg peaks (vertical bars). The two panels show 
the XRD pattern obtained (a) before and (b) after annealing 
the sample. 
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FIG. 15: (Color online) Room temperature powder XRD pat- 
tern (red circles) of LaNi2(Po.75Geo.25)2, Rietfeld refinement 
fit (solid black line), difference profile (lower solid blue line), 
and positions of Bragg peaks (vertical bars). 
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FIG. 16: (Color online) Room temperature powder XRD pat- 
tern (red circles) of LaNi2(Po.5oGeo.5o)2, multi-phase Rietfeld 
refinement fit (solid black line), difference profile (lower solid 
blue line), and positions of Bragg peaks (vertical bars; upper: 
LaNi2(Po.5oGeo.5o)2, lower: NisGes). The refinement reveals 
that the phase composition is 99.6 mol% LaNi2(Po.5oGeo.5o)2 
and 0.4 mol% NisGes. 
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FIG. 17: (Color online) Room temperature powder XRD pat- 
tern (red circles) of LaNi2(Po. 25060.75)2, Rietfeld refinement 
fit (solid black line), difference profile (lower solid blue line), 
and positions of Bragg peaks (vertical bars). 
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FIG. 18; (Color online) Room temperature powder XRD 
pattern (red circles) of LaNi2Ge2, multi-pfiase Rietfeld re- 
finement fit (solid black line), difference profile (lower solid 
blue line), and positions of Bragg peaks (vertical bars; up- 
per: LaNi2Ge2, middle: La203, lower: Nii9Gei2). The re- 
finement indicates that the phase composition of the sam- 
ple is 93.9 mol% LaNi2Ge2, 5.7 mol% La203, and 0.4 mol% 
Nii9Gei2. 
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FIG. 19: (Color online) Magnetization M versus applied field 
H isotherms (symbols) at the listed temperatures for LaNi2P2 
(a) before and (b) after annealing. The solid curves of the 
corresponding colors are fits by Eq. (19). 



Appendix C: Magnetization versus Field Isotherms 

In this Appendix the M{H) isotherms in Figs. 19-21 
and the dependence of the ferromagnetic impurity satu- 
ration magnetization on temperature in Fig. 22 that are 
discussed in Sec. IV C are shown. 
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FIG. 20: (Color online) Magnetization M versus applied mag- 
netic field H isotherms (symbols) at the listed temperatures 
for LaNi2(Gei_2:Pa;)2 samples with (a) x = 0.25, (b) x = 0.50, 
and (c) X — 0.75. In each figure, the solid curves of the cor- 
responding colors are fits by Eq. (19). 
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FIG. 22: (Color online) Saturation magnetization of ferro- 
magnetic impurities Mo in the LaNi2(Gei-2:Pa;)2 samples ver- 
sus temperature T . The data at 350 K are extrapolated. The 
data for x = 0.5 suggest the presence of a ferromagnetic im- 
purity phase with a Curie temperature of « 100 K. The solid 
lines are guides to the eye. 
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FIG. 21: (Color online) Magnetization M versus applied 
field H isotherms (symbols) at the listed temperatures for 
LaNi2Ge2 (a) using a SQUID magnetometer and (b) using a 
VSM magnetometer. The solid curves of the corresponding 
colors are fits by Eq. (19). 
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